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0> ■ Abstract 

f~^ \ We analyze the known moving boundary effects in plasma and briefly review a 

mechanism for particle acceleration through a plasma front. We then show how this 

mechanism may aff^ect particles crossing the boundary of an expanding quark gluon 

i-Q \ plasma. 

P-i! 

Oh: 

> 

X 



1 Introduction 



Heavy ion collisions have been studied all over the years, and their current interest lies 
mainly on the search for the Quark-Gluon Plasma (QGP) state |^. Information on the 
eventual transition of the hot nuclear medium to the plasma state can be obtained from 
the emission of energetic particles (electrons and muons) resulting from the heavy ion 
colhsion [0. 

Usually it is assumed that particles coming out of the plasma region have their energy 
equal to the one they had inside the QGP blob. However, due to the collective interactions 
with the surrounding plasma medium, the final particle energy can significantly change 
when they move across the boundary of the expanding plasma. This energy shift is 
due, not only to the difference in the value of the effective potential created by all the 
surrounding particles on each side of the boundary, but also to the velocity of the moving 
boundary. 

This collective effect of the moving background medium has been explored by several 
authors in a different context, and constitutes the basis for plasma acceleration. In this 
respect, photon acceleration by relativistic plasma perturbations is a well defined concept 
0, which has been tested experimentally with success [Q. A similar problem, suggested 
by Bethe 0, involves the neutrino plasma interactions, where the electromagnetic force is 
replaced by the weak force and can eventually be relevant to dense plasmas in supernovae 
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From the point of view of the QGP which might have, or has already, been produced 
in high energy heavy ion collisions, this type of collective processes can have a direct 
consequence, for instance, on the behaviour of dileptons produced in the plasma region. 
For them, QGP is a normal high temperature electromagnetic plasma and therefore one 
can approach the process with techniques previously developed for this kind of medium. 



2 Energy states across a plasma boundary 



Let us start with the dispersion relation of a particle moving in a QGP (in units c = h = 1): 

^E-Vf=p^ + m^ (1) 

where p is the particle momentum, E its energy, m its mass and V is the effective potential 
created by the interactions with the surrounding particles of the plasma medium. From 
here we can obtain the wave equation for the particle normalized wave function: 

(v^_„._g,.2.^|_^V (2) 

We shall be interested in a plasma density perturbation moving with a constant ve- 
locity Vg. For simplicity, we assume that the perturbation is uniform in the plane perpen- 
dicular to Vg. We can then choose the x coordinate parallel to that velocity. Using the 
d'Alembert transformations to go to the moving reference frame: 

$, = X — Vet , rj = X + Vgt (3) 

and performing a Fourier transformation in 7], defined as: 

<l>{^,v) = jMOe''''^ (4) 

we obtain: 



(1 - ^') (^ - <l") + 2^?(1 + ^s)^ - ^' 



0.(0 



= 2zqVgV{0 (tq - ^] MO + V\OMO (5) 

If we consider that we are in an equilibrium QGP phase, then we must accept that 
the lifetime of the system is larger than the relaxation times associated with the processes 



occurring in the medium. In particular, we must assume that the potential V{^) is a 
slowly varying function. If this is true, we can use the WKB approximation and set: 



MO = 0goexp iij piOd^j (6) 

^From this we can get the local dispersion relation: 



(1 - v^,) [p\0 + q'] + 2(1 + v^,)q + m' = 

= 2vsV{0[q-pm + V\0 (7) 

We now take the usual definitions: 



k^-i^ , £ = ^ = 4 (8) 

We can then establish the relations between the local frequency (or energy) uj{$,), the 
local momentum k{C,), the local variable p(^) and q: 

2p{o = m + ^ (9) 



2g = kiO - ^ (10) 

This allows us to write the dispersion relation (^ in the form: 



P{0 +m'- co^O = -2coiOViO + V\0 (11) 

And, solving for uj{^), we get: 



u;{0 = V{0±Jk^{0 + m' (12) 



Now, by definition, g is a constant parameter, which for a given quantum state is 
invariant inside and outside the plasma region. Thus, we can quite easily relate the initial 
values of the energy and momentum, tUj = uj^C^i) and ki = k{^i), for a particle created 
inside the plasma, with the final values ujf and kf, observed outside: 

2qvs = kiVs -uJi = kfVs - ujf (13) 

The energy shift of the particle will then be: 

Acu = ojf — uji = {kf — ki)vs (14) 

But the initial and final momenta, ki and kf can also be expressed in terms of the 
energies, according to equation (pT]). Assuming that the effective potential is equal to 
zero outside the plasma, V/ = 0, because there are no interacting particles, we get: 



ki = uji\/{l - Vi/uiY - IV? /uj?, 

kf = ujf^l - m^/uj} (15) 

Replacing this in equation ([TsD , we obtain: 



^f = ^^ /. ,,2 (^6) 

1 — VsJl — Vfl'^ jUJf 

In order to get an explicit result for the final energy, let us assume that the rest mass 
of the particle m is negligible, which is certainly true for the electrons. The results is 
then: 



UJf ~ UJi 



\-v,{\-Viluj,) 



Which means that the energy shift (|14|) is simply given by: 



(17) 



Auj ~ uji '"' (18) 

This means that, when Vg tends to 1, the energy shift can be considerable, even if 
the initial effective potential inside the plasma is a small correction to the energy of the 
ejected particle, Vi ^ cjj. 

Notice that if Vg < 0, that is, if we have an imploding surface, leptons traversing it 
will lose energy. For f s ~ 1 leptons can change their energy by as much as half their value 
inside the plasma region. 



3 Particle dynamics 



A few remarks should be made concerning the above result. First, equations (|l^ - |T8D 
are only valid for Vf = 0, which means that the particle has to leave the plasma region 
before it spreads away. Second, the exactly resonant case, corresponding to f^ = 1 is 
not determined by these equations because in this case the particle will never cross the 
boundary and, as a result, its energy remains constant: Auj = 0. Finally, the above 
model of a constant velocity Vs can only be relevant if the changes in the plasma main 
parameters, such as density and temperature, are negligible during the interaction time 
with the moving boundary. 

These, and other, aspects of the problem can be elucidated by using the classic equa- 
tions of motion for the emerging particle. They can be used in the same spirit of the 
WKB approximation. These equations can be stated as: 



df duj dk duj 

'd^^'dk ' ~di^"df ^^^^ 

Here, r, k and u are the mean values of the position, momentum and energy of the 
particle. In a one dimensional model, we can write: 



dx 


dk 


= 




2k 


dt 


Vk 


2 _|_ ^2 


dk 
dt 


= — 


did 
dx 


= 


dv 

dx 



(20) 



(21) 

Notice that the energy ui will change when the particle moves across the boundary, 
according to: 



di^^du^^dV_ 

dt dt dt ^ ' 



For a boundary moving without deformation with velocity Vg, we can use V{x,t) 
V{x — Vst), which leads to: 



^ = -1^ (23) 

dx V. dt 



Comparing with equations ( pl| , p^), we conclude again that {u — kvs) is a constant of 
motion, as considered before. 

Let us consider the case where V{x,t) 7^ V{x — Vst), and this constant of motion is 
not valid. As a first example we can use: 

Vix, t) = Voit)^ {1 + tanh [hix - Xo(t)]} (24) 

where k^ defines the slope of the boundary, and Xo(t) = xq + Vst is its mid position. Here, 
we can also define the amplitude of the potential Vo(t) as nearly proportional to the mean 
plasma density. Actually, there is also a small dependence on the plasma temperature, 
associated with Debye screening. If the plasma is initially inside a sphere of radius vq, 
and it is expanding along the axis x, we can then approximately write: 

Voit) = V. \ (25) 

1 + {Vs/ro)t 

The energy shift of the particle moving across the boundary along the x axis can be 
determined by: 

uj{x, t) = V{x, t) + yP(t)Tm2 (26) 

where the particle postion x{t) and momentum k{t) are determined by the above equations 
of motion. 

A more realistic example is provided by the expanding fireball model P] for which we 
assume a cylindrical expanding plasma region such that the volume evolves like 

Voit) = 277 |xo + Vst + -a,.t'| (ro + -a^t^ (27) 
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where the parameters are fixed using fiow data. Typical values are tq ~ 4.6 fm an a 
total freezout time of tj ~ 10 — 12 fm in the longitudinal x direction and 3 — 4 f m along 
the perpendicular direction, which implies a final velocity of the front along the x axis, 
Vx ~ 0.75c and v± ~ 0.55c. In this case, for two particles with the same energy inside the 
plasma Ui, but crossing different plasma boundaries we have the relation 



c^/ii l-Vs±' l-Vs\\{l-Vi/uJi) 



{2i 



In order to illustrate the relevance of the collective plasma effects and to point out the 
main aspects of the problem of determining the initial energy of the observed particles, 
we present numerical examples resulting from the integration of the equations of motion. 
In figure 1, we show the energy shift for two different values of the front velocity Vg-, and 
a constant plasma potential. We notice that we can obtain a non-negligible energy shift 
AciJ ^ V , which increases when Vg approaches one. This illustrates the importance of 
the collective effects associated with a rapidly moving potential on the particle dynamics. 
In figure 2, we show that such energy shifts can be reduced if the maximum potential 
amplitude decreases during the interaction time with the boundary. Even so, there is still 
a large range of parameters for which the energy shifts can be significant. 
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4 Conclusions 



We have shown that the energy of a particle, for instance, an electron or a muon, measured 
during high energy ion collisions and emerging from the dense nuclear matter gas or quark 
gluon plasma can be significantly different than its initial energy. This results from the 
influence of the effective potential acting on that particle as a consequence of its interaction 
with the background plasma particles. Even if this effective potential remains negligible 
with respect to the initial particle energy, its influence can be important due to the 
fact that it is changing on very fast space and time scales. In particular, the potential 
time variation create a force which attains its maximum value at the expanding plasma 
boundary. Even if this force is weak, it can act for a long time if the particle moves along 
with the boundary. The results are valid for electromagnetic and for strong interactions 
with the background particles. They suggets that the mean field acceleration processes 
already known in other other areas of plasma physics can also be relevant to high energy 
physics. 

This work shows that the resulting energy shifts are very sensitive to the space and 
time variations of the expanding plasma, and in particular to the velocity of its boundary. 

In this paper we have used the WKB approximation, which is not valid if the back- 
ground plasma potential changes on scales comparable to the wavelength of the emerging 
particles. However, it can easily be realized that a full quantum model will lead to the 
same energy shifts. Typical quantum processes not described by the present model, such 
as reflection of particles from the boundary and coupling between different helicity states 
will be discussed in a future work. 
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FIGURE CAPTIONS 

Figure 1 - Particle energy as a function of time, for rn^ = 0.1, for a tangent hyperbolic 
front with ks = 1, and a constant potential Vi = 1. The front velocities are: (a) Vs = 0.98, 
and (b) Vs = 0.99. 

Figure 2 - The same as in Figure 1, but for a time dependent potential, as defined by 
equation (pSf), with Vs/tq = 0.1. 
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